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Abstract. Understanding species responses to climate change is crucial for the mitigation of its effects. Few studies, however, 
have examined how climate change impacts timing in reptile life cycles, or how it may interact with other life history traits. 
Here, we explore associations between climatic variation, reptile phenology, and juvenile growth. We used the chameleon spe- 
cies Furcifer labordi as a biological model given its unusual annual life history: embryonic diapause during austral winter is 
followed by synchronous hatching and an extremely short adult lifespan. We sampled chameleons (n = 674) and climate data 
across three field seasons between 2013-2018 in Kirindy Forest (CNFEREF), western Madagascar. Evidence of temperature and 
precipitation shifts have already been reported for this region. We show climate-dependent variation in the timing of hatching: 
it always occurred after the first rains and its onset varied by up to 32 days among years (equivalent to 12-17% of the total length 
of non-embryonic lifespan). Moreover, delayed rainy season onset resulted in a shortening of the period between the first rains 
and first hatching. Significant differences in hatching time and temperature also resulted in variation in juvenile growth, spe- 
cifically faster post-hatching growth in the late-hatching cohort. Therefore, climatic variation was associated with phenological 
and juvenile growth plasticity in the life cycle of E labordi. Implications for conservation are worrying, as complete populations 
could be lost if a single cohort of this annual species cannot effectively respond to increasing climatic fluctuation. 
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Introduction 


Among vertebrates, poikilothermic taxa such as reptiles are 
hypothesized to be particularly dependent on environmen- 
tal conditions for growth and reproduction. Together with 
their often-limited dispersal abilities, this constraint makes 
reptiles particularly vulnerable to rapid environmental shifts, 
and they may suffer higher extinction rates than other taxa 
due to climatic change (SINERVO et al. 2010, Huey et al. 
2012). In various reptiles, less dramatic effects of tempera- 
ture changes have been observed for hatchling phenotypes 
(e.g., BOOTH 2006, SINGH et al. 2020), hatchling behaviour 
(DAYANANDA et al. 2017, SIVITER et al. 2017), hatching suc- 
cess (Diaz-PANIAGUA et al. 2006), post-hatching survival 
rate (ANDREWS et al. 2000), and, in some species, offspring 


sex (e.g., JANZEN & PAUKSTIS 1991, SCHWANZ et al. 2010, BAL- 
LEN et al. 2016). 

In addition, range shifts due to climatic variation have been 
repeatedly documented in reptiles, including lizards (HouNI- 
ET et al. 2009, Ce1A-HassE et al. 2014), snakes (LOURENGO- 
DE-Mo rags et al. 2019), and chelonians (IHLow et al. 2012). 
However, many populations now exist in isolated habitats 
without access to other potentially suitable habitats, making 
range shifts an unlikely response to less favourable environ- 
mental conditions. It is therefore important to better under- 
stand in situ responses to climatic variability across a wide 
range of taxa. Such information is interesting for understand- 
ing plasticity in life history traits and crucial for identifying 
species threatened by climate change for ecological forecast- 
ing and conservation planning (AKÇAKAYA et al. 2006). 
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Among in situ responses, modification of the timing of 
events in a life history can occur in response to abiotic fac- 
tors exhibiting variation over time (egs temperature, pre- 
cipitation) (BELLARD et al. 2012). Such phenological shifts 
are well documented in many taxa (e.g., STEFANESCU et al. 
2003, RATHER et al. 2018, STEPANIAN & WAINWRIGHT 2018, 
STAUDINGER et al. 2019). Studies focusing on reptiles’ phe- 
nological responses to climate change are scarce, but few 
publications (e.g., MORENO-RuEDA et al. 2009, CLARKE & 
ZANI 2012, LJUNGSTRÖM et al. 2015), focus on the shifts of 
life history traits such as growth pattern caused by higher 
ambient temperatures. 

Here, we study variation in the phenology and growth 
in annual cohorts of Labord’s chameleon (Furcifer labordi) 
across years with significant climatic variation. This species 
is characterized by an extreme life history among terrestri- 
al vertebrates as it is thought to have the shortest lifespan 
(4-9 months) among tetrapods as part of an unusual an- 
nual life cycle (KARSTEN et al. 2008, ECKHARDT et al. 2019). 
These chameleons are almost exclusively represented in the 
egg stage (embryonic diapause) during the austral win- 
ter followed by synchronous hatching, extremely rapid 
growth, and reproduction before the onset of the next dry 
season (KARSTEN et al. 2008). As an entire population con- 
sists of a single, synchronous age cohort, F labordi provides 
a unique opportunity to monitor the effects of climate fluc- 
tuations on phenology and growth for entire generations 
of a tetrapod species. Further, due to the extremely com- 
pressed lifespan, perturbations to timing could substan- 
tially reduce the time window available for growth and re- 
production, raising questions about which components of 
the life history are more or less plastic. 

Furcifer labordi is restricted to dry forest patches in 
western and southwestern Madagascar and is classified as 
Vulnerable in the IUCN Red List (IUCN 2021). Two pop- 
ulations have been studied: a northern population in Ki- 
rindy Forest, within the protected area of Menabe Anti- 
mena (ECKHARDT et al. 2017, 2019), and a southern pop- 
ulation in the protected area of Ranobe (KARSTEN et al. 
2008). Here, we focus on the northern population from 
Kirindy, where the abundance of F. labordi is still high 
(RANDRIANANTOANDRO et al. 2010) and the effects of cli- 
mate change have already been observed. GoopMan et al. 
(2018) analysed climatic fluctuations in the protected area 
of Menabe Antimena over 30 years (1985-2014). An annual 
increase of approximately 0.5% in precipitation has been 
observed and the rainy season shifted on average 10 days 
later in more recent years. In addition, the average mini- 
mum and maximum daily temperatures increased by 1.4°C 
and 0.7°C, respectively. 

In this study, we estimate the timing of emergence and 
the growth rates of hatchlings and juveniles for three gen- 
erations of F. labordi, spanning years with significant cli- 
matic variation. First, we hypothesized that a delayed rainy 
season onset would be associated with a delayed hatching 
date. Second, we hypothesized that juvenile daily growth 
rate increases in late-hatching cohorts to compensate for 
this delay. 
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Materials and methods 
Study site and species 


The present study was carried out in the dry deciduous Kir- 
indy Forest (CNFEREF) (44°28’-44°46’ E, 20°03’-20°10’ S, 
30-60 m a.s.l.), in the Harmonious Protected Landscape 
of Menabe Antimena, western Madagascar. Although orig- 
inally covering 12,500 ha, the integrity of Kirindy Forest 
is currently threatened due to habitat loss and fragmenta- 
tion (GOODMAN et al. 2018). The overall rate of deforesta- 
tion within the protected area was 0.67% per year with a 
maximum of 2.55% between 2008 and 2010 (ZINNER et al. 
2014), presumably associated with the 2009 political crisis 
in Madagascar. 

Regarding local climatic conditions, the study area lies 
in a transitional zone. The seasonally dry climate typical 
of the west coast of Madagascar dominates in the north- 
ern part of the protected area. In contrast, the sub-arid cli- 
mate typical of southwest Madagascar dominates in the 
southern portion (GOODMAN et al. 2018). Rainfall averages 
954mm per year (1981-2017), of which 96% falls between 
November and April. Temperature ranges from a mean 
daily minimum of 20.0°C to a mean daily maximum of 
30.3°C. June to August is the coolest period (min, 12.4°C) 
and September to November the hottest (max. 33.4°C). Ki- 
rindy Forest is characterized by intra-annual variation (i.e., 
extreme seasonality) and unpredictable inter-annual var- 
iation in climate (KAPPELER & FICHTEL 2012), providing 
an opportunity to study the effects of substantial climatic 
changes over short time scales. 


Vegetation cover measurement: 
Normalized Difference Vegetation Index (NDVI) 


NDVI, calculated by the reflection rate of red light (RED) 
and near infrared light (NIR) [NDVI=(NIR-RED)/ 
(NIR+RED)], is strongly correlated with aboveground net 
primary productivity (MyNENI et al. 1995; see also http:// 
daacmodis.ornl.gov). Here, NIR and RED are the rates of 
near infrared light and red light that are reflected by the 
vegetation and recorded by the satellite sensor. The for- 
mula results from the fact that chlorophyll absorbs RED, 
whereas NIR is re-emitted. To compare seasons, we meas- 
ured the proportion of vegetation cover in the middle of 
the dry season (Jul.-Aug.), the late dry season (Sep.—Oct.), 
and the wet season that covers the hatching time and pe- 
riod of juvenile growth (Nov.—Dec.) in 2013, 2015 and 2017. 
For analysis of NDVI, we selected a square of 3km” around 
the centre of our study site. 


Climate data analysis 


Temperature and precipitation data were collected from 
2013 to 2018 automatically each day, every 10 minutes, with 
a weather station (Lambrecht meteo GmbH, Göttingen, 
Germany) at the research camp in Kirindy Forest. We ex- 
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amined temperature and precipitation fluctuations during 
the field seasons using Kruskal-Wallis test, and visualized 
the variation via heatmaps. 


Chameleon sampling 


Chameleon sampling took place over three field seasons 
that cover the hatchling and juvenile growth phases in 
the E labordi life cycle: 18 November 2013 to 14 February 
2014; 12 October 2015 to 17 December 2015; and 11 Novem- 
ber 2017 to 20 December 2017, 16 January 2018 to 14 Feb- 
ruary 2018. Chameleons were searched using a headlamp 
and a Maglite torch at night as they perch on distal ends 
of vegetation and display a conspicuous pale or vivid col- 
our when lit (JENKINS et al. 1999). Two different transects 
were surveyed alternately between 7 pm and midnight. In- 
dividuals were captured by hand or removed from higher 
branches using an extendable pole or stick to which they 
were encouraged to grip. We scheduled a recurring order 
of 10 nights of sampling followed by a break of four nights 
within each field season. 

For each individual, capture location was marked with 
an individually numbered coloured flag and coordinates 
were collected using a GPS. Age class (hatchling, juvenile, 


or adult) was recorded using the following features: (i) in- 
dividuals with a distinct small opening in the abdominal 
wall were classified as hatchlings, (ii) individuals lacking 
the abdominal opening but not yet sexually mature were 
classified as juveniles, and (iii) sexually mature adults were 
recognized by the appearance of a colourful mating pat- 
tern in females and the presence of a distinct hemipenis 
bag (bulge at the base of the tail) and a hard, ossified ros- 
tral appendage in males. The rostral appendage is already 
present in the juvenile stage but is not fully developed until 
the adult stage (Fig. 1). The abdominal opening in hatch- 
lings is caused by a former connection with the yolk sack 
and is visible for about one week after hatching. Sex was 
determined in juveniles based on the presence/absence 
of the rostral appendage, and in adults based on the pres- 
ence/absence of the hemipenis bag and rostral appendage. 
Snout-vent length (SVL) was measured with a calliper to 
the nearest 0.1 mm. All individuals were released at the ex- 
act point of capture following manipulation. 


Hatching date and climate 


We estimated the date of hatching of F. labordi for each 
year from the first date hatchlings were observed during 


Figure 1. Individuals of Furcifer labordi in life: (A) hatchling climbing over a juvenile male that is several weeks older; (B) juvenile 
female short before change to adult coloration; (C) adult male and female during mating; (D) adult female displaying stress coloration. 
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the nocturnal surveys and from estimates of growth rates 
for hatchlings and juveniles (described below). To account 
for the possibility that unobserved hatching occurred pri- 
or to the first observation of hatchlings, we use growth 
models to determine the probable number of days since 
hatching for juveniles of a given size. To examine the re- 
lationship between timing of hatching and rainfall in this 
population, we compared the estimated date of first hatch- 
ing with the start date of the rainy season (as defined by 
the first date with measurable rainfall) for 2013, 2015, and 
2017. 


Analysis of growth patterns 


Due to fluctuations in body mass with variation in stom- 
ach contents, fat storage, and reproductive status in lizards, 
length is more appropriate to measure growth than weight 
(ANDREWS 1982, HAENEL & JOHN-ALDER 2002). Therefore, 
our analysis of growth focuses on SVL. For each field sea- 
son, we defined the date of observation of the first hatch- 
ling as day zero. After checking the assumptions of para- 
metric tests, growth was estimated by fitting log(SVL) (re- 
sponse variable) on days since day zero (predictor variable) 
for each year’s cohort of hatchling and juvenile chamele- 
ons. Data analyses and data visualization were carried out 
using the software R 4.0.5 (R Core-Team 2021) and the sig- 
nificance level was set at P < 0.05. 


Results 
Vegetation cover measurement: 
Normalized Difference Vegetation Index (NDVI) 


Comparisons across years for a given season revealed that 
2013 and 2017 were consistently characterized by less veg- 
etation cover than 2015, with 2013 showing the lowest val- 
ues in general (Table 1). The higher values of vegetation 
cover in 2015 probably result from heavier rainfall in the 
previous rainy season. Reviewing average vegetation cover 
from 2008 until 2017, we found that the vegetation cover 
calculated for 2013 and 2017 was around average in the dry 
season and at the end of the dry season, whereas 2015 was 
above average. Concerning the vegetation cover during the 
beginning of the rainy season (main time of hatching), 2013 
and 2017 were below average, whereas 2015 was above av- 
erage. 


Climatic variation 


The average daily temperature (Kruskal-Wallis: x? = 
282.04, df = 4, P < 0.001) and precipitation (Kruskal-Wal- 
lis test: x? = 117.91, df = 4, P < 0.001) differed significantly 
during the study period (Fig. 2). The lowest monthly mean 
minimum temperature (23.7°C) was recorded in December 
2015, whereas the highest monthly mean maximum tem- 
perature (30.28°C) was recorded in January 2018. 
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Table 1. Averaged NDVI for the onset of the dry season (Jul.- 
Aug.), end of the dry season (Sep.—Oct.), and the rainy season 
covering hatching and juvenile growth period (Nov.—Dec.). 


Year Jul./Aug. Sep./Oct. Nov./Dec. 

2013 0.6153+0.04 0.5888+0.03 0.5945+0.05 
2015 0.7197+0.03 0.641+0.01 0.6864+0.09 
2017 0.6633+0.03 0.5796+0.02 0.6213+0.09 
2008-2017 — 0.6309+0.05 0.5677+0.05 0.6415+0.125 


Chameleon sampling 


During the field surveys conducted for this study, we sam- 
pled 674 individuals of Furcifer labordi. Of these, 73 (11%) 
were hatchlings, 421 (62%) were juveniles, and 181 (27%) 
were adults (Table 2). Hatchlings were observed from mid- 
October to December, except one hatchling found in Jan- 
uary 2018, most likely an outlier. Juveniles were mainly 
found from November to January. Adults were observed 
from January onwards, apart from one adult female found 
in October 2015, most likely a survivor from the previ- 
ous active season (Fig. 2). Indeed, although the adults of 
F. labordi generally die at the onset of the relatively cold 
dry season, a small number of females might survive until 
a second breeding season when the rainy season is unu- 
sually long (ECKHARDT et al. 2017). In terms of sex ratio, 
we found a male-bias in adults for the 2013-2014 cohort 
(Chi-squared: y” = 5.26, df = 1, P < 0.05). Males are sig- 
nificantly larger than females in both juveniles (Wilcoxon 
test: W = 14076, P < 0.05) and adults (W = 69, P < 0.05) 
(Table 2). 


Hatching time variation 


The first hatchlings were observed on 20 November in 
2013, 17 October in 2015, and 18 November in 2017 (Fig. 2). 
However, in 2013, juveniles were observed on the same date 
as the first hatchlings. As an almost exclusively annual spe- 
cies (with the rare exceptions noted above), populations of 
F. labordi consist of single, yearly cohorts such that gen- 
erations do not overlap. This means the juveniles seen on 
20 November 2013 likely indicate the existence of earlier 
hatching that occurred prior to the start of sampling in 
2013 (on 18 November). Using the estimated average daily 
growth derived from growth models (Fig. 3) and the SVLs 
of the first juveniles observed, we estimate that hatching in 
2013 likely occurred around 27 October. 

In 2015 and 2017, the first observations of hatchlings 
were preceded by more than 5 days of careful searching 
that resulted in no observed hatchlings. The first juveniles 
were observed 18 to 31 days after the observation date of the 
first hatchlings. From growth models, the estimated date of 
hatching is 13 October to 22 October, and 1 November to 18 
November, respectively, for 2015 and 2017. 
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Table 2: Number and mean snout-vent length (SVL) in mm of detected Furcifer labordi individuals per age class and sex in 2013, 


2015, and 2017. Sex not determined for 8 juveniles (NA). 


Juveniles Adults 
Year n  Hatchlings F M NA F M 
2013-2014 362 14 119 144 0 30 55 
2015 142 46 44 43 8 1 0 
2017-2018 170 13 35 28 0 44 50 
SVL (meant SD) 674 26.642.3  44.6412.7 51.94+13.6 NA  73.5+3.3 94.0+9.3 


Hatching time and climate 


In order to investigate the relationship between hatching 
time and climate, the date of hatching was mapped onto 
the onset of the rainy season (Fig. 2). Hatching always oc- 
curred after the first rains. The rainy season started on 
2 October in 2013 (precipitation = 25.7mm), and the es- 
timated date of hatching (27 October 2013) occurred 25 
days after first rains. The interval increased in 2015, with 
40 days after the first rains, which started on 8 Septem- 
ber (precipitation = 6.9 mm). However, in 2017, hatchlings 
were observed just one day after the first rain on 17 No- 
vember (precipitation = 2.9 mm). As with the hatching 


later the start of the rainy season, the shorter is the pe- 
riod between first rains and first hatching. Additionally, 
the amount of first precipitation decreased over the stud- 
ied years. 


Growth rate 


Hatchling and juvenile growth followed a linear regression 
model (R? = 0.66, P < 0.001), and the regression equations 
for each cohort are shown in Figure 3. Individual growth 
rate was significantly higher for the cohorts from 2013- 
2014 and 2017-2018, compared to the cohort from 2015 


time, the start of the rainy season varied across years. The (ANOVA: F ay = 75 P< 0.001). 
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Discussion 


The results of our study showed climate-dependent varia- 
tion in the timing of hatching of F labordi, in which a de- 
layed onset of the rainy season resulted in a shortening of 
the period between first rains and first hatchings. Addi- 
tionally, significant differences in hatching time and am- 
bient temperature caused differences in juvenile growth, 
in which faster post-hatching growth occurred in the late- 
hatching cohort. 

Previous studies of the southern population of E labordi 
in Ranobe showed a synchronous hatching at the onset of 
the rainy season (KARSTEN et al. 2008). Our results from 
the northern population in Kirindy showed that hatching 
occurred after the first rains and hatchlings can generally 
be detected for over a month. A single hatchling was also 
found at the end of December in 2013 and around mid- 
January in 2018. We therefore observed an extended hatch- 
ing period in Kirindy. As the hot rainy season is longer in 
Kirindy compared to Ranobe (GoopMan et al. 2018), food 
resources are more abundant and available for a longer pe- 
riod. 

We also observed variation in the hatching time of 
F. labordi that follows a pattern similar to the rainfall 
schedule. The hot rainy season started earliest in 2015, fol- 
lowed by 2013 and 2017, and the earliest observed hatch- 
ing of F. labordi followed the same order. Based on growth 
models, estimated hatching time only differed by 10 days 
between 2013 and 2015. Moreover, we found that first ju- 
veniles appeared around the same period. However, over 
a month separated the hatching time in 2017 compared 
to 2015. Juveniles were also observed approximately three 
weeks apart, actually indicating delayed hatching in 2017. 
The relatively long time span between the detection of the 


y = 3.379 + 0.0121 x R? = 0.63 
y = 2.941 + 0.0148 x R? = 0.55 
y= 3.330 + 0.0144 x R? = 0.83 
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Figure 3. Graph showing linear regression model of hatchling and 
juvenile growth of Furcifer labordi on time in Kirindy Forest. De- 
picted is log-transformed snout-vent length (log(SVL)) against 
days post hatching for the three cohorts (2013-2014, 2015 and 
2017-2018). 
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first hatchlings and juveniles might be due to the fact that 
hatchlings from one clutch emerge synchronously and 
sleep in the proximity of the nest in the first night(s) and 
are therefore easier to detect than those that already dis- 
persed. Additionally, at the beginning of the rainy season 
the probability to detect small juveniles is rapidly decreas- 
ing due to increasing foliage growth. 

Delayed hatching can occur when environmental con- 
ditions are not suitable (Doopy, 2011). A delayed rainy 
season accompanied by a decreased amount of rainfall 
in 2017 have probably resulted in the later emergence of 
hatchlings. The eggs of chameleons are sensitive to eco- 
logical conditions during incubation, especially variation 
in temperature and moisture (DÍAZ-PANIAGUA et al. 2002, 
Diaz-PANIAGUA & CUADRADO 2003). In fact, their eggs are 
flexible-shelled as they contain a minimal calcareous layer 
(PACKARD et al. 1982), making them more permeable and 
more vulnerable to hydric conditions during incubation 
(e.g., PACKARD 1991, 1999). As the eggs are buried under- 
ground and rainwater filters through the soil to reach the 
nest, the soil should be moist enough for hatchlings to be 
able to crawl and find a way through to reach the surface 
(MEAsEY et al. 2014). 

Henceforth, the observed variation in the onset of the 
rainy season induces a temporal shift in the hatching time 
of F. labordi. Nevertheless, as temperature also plays a role 
in the activation of embryonic development and ending of 
diapause (ANDREWS et al. 2008), there is also a need to in- 
vestigate the relationship between timing of hatching and 
temperature fluctuation. The observed changes in the phe- 
nology of F. labordi may indicate that climate change is af- 
fecting living systems in negative ways (e.g., PARMESAN & 
YOHE 2003). At the population level, the observed short- 
ened interval between the onset of the rainy season and the 
hatching date from 25 days in 2013 to one day in 2017 sug- 
gests a possible risk of future mismatches between the two 
events, resulting in decreased hatching success. 

Climate change might also have additional effects on 
this species. For instance, within the population in Ki- 
rindy, male adult body size differed significantly within 
years (ECKHARDT et al. 2017). This observation was linked 
to inter-annual climate fluctuations, suggesting a relatively 
fast adaptation to shorter reproductive seasons concern- 
ing adult body size. This assumption is underpinned by the 
observation that males in Kirindy were significantly larg- 
er compared to those in the population of Ranobe (Ecx- 
HARDT et al. 2017). Therefore, the reduction of adult body 
size might be an answer to climate change in this species as 
well. Moreover, BALLEN et al. (2016) observed that the sex 
of hatchlings in the Veiled chameleon Chameleo calyptra- 
tus is affected by the interaction of egg mass and tempera- 
ture. This phenomenon might be more widespread within 
chameleons and climate change might additionally lead to 
a sex bias in F labordi. ECKHARDT et al. (2019) observed a 
strong sexual bias towards males in juvenile F nicosiai that 
might be ascribed to partial temperature sex determina- 
tion in this species. At the community level, three sympat- 
ric Furcifer species are observed in the Kirindy Forest and 
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they probably partition the resources by the temporal shift 
in species-occurrence of hatching, with F labordi hatch- 
ing the earliest (ECKHARDT et al. 2019). Therefore, the ob- 
served delayed hatching of F labordi may increase compe- 
tition for food and space resources between the three sym- 
patric species. However, further analyses are required to 
clearly document the impacts of climate change on species 
interactions (e.g., FORCADA et al. 2006, YANG & RUDOLF 
2010). These are preliminary results, but long-term studies 
with experimental demonstration are necessary to obtain 
results that are more conclusive. Moreover, data from ad- 
ditional generations are needed to determine if these pat- 
terns hold across a larger sample of years. Furthermore, the 
shift in the onset of the rainy season and precipitation pat- 
terns that we observed might be highly linked to the In- 
dian Ocean Dipole (IOD) that is responsible for irregular 
oscillation of the sea surface temperatures, but long-term 
recordings are obligatory to investigate a clear correlation. 
The IOD has only been recognized about two decades ago 
(WEBSTER et al. 1999), making it much less studied than 
other climate patterns. However, increased frequency of 
extreme IOD events due to greenhouse worming are pre- 
dicted in the study of Car et al. (2014), leading to extreme 
climate and weather events in regions that are affected by 
the IOD, such as Madagascar. 

In general, short-lived species invest more energy into 
growth and reproduction at the expense of health mainte- 
nance (e.g., SMITH & FRENCH 2017). We confirmed the pre- 
dicted increased investment in individual growth to com- 
pensate for delayed hatching in this study. Several other 
species also exhibit faster post-hatch growth in late hatch- 
ers (e.g., BENOWITZ-FREDERICKS & KITAYSKY 2005, NEW- 
TON 2008, WALLUS & SIMON 2008, BALDASSARRE 2014). A 
faster growth is necessary for these hatchlings to reach the 
stage of sexual maturation quickly. As observed through- 
out the field seasons, gravid females appeared around the 
same period despite the observed variation in hatching 
time. Thus, juvenile growth provides a mechanism to flex- 
ibly adjust the life history schedule of late-hatched cohorts, 
but the clock they use to calibrate their short life cycle re- 
mains elusive. In addition, the main egg-laying period 
likely remains constant, whereas elongated rainy seasons 
can lead to extended egg-laying periods (ECKHARDT et al. 
2017). 

Concerning the average vegetation cover within the 
growth period of F. labordi, 2015 showed the highest val- 
ues, indicating that more water and insect food were avail- 
able for the growing juveniles. However, in this season, ju- 
veniles grew more slowly compared to the other two sea- 
sons. The daily temperature and precipitation varied signif- 
icantly during the study period. Variation in growth is as- 
sociated with climatic conditions, especially temperature, 
and the effects are particularly pronounced in reptiles. For 
example, juvenile snapping turtles (Chelydra serpentina) 
grow slowly in a cool (19°C) versus a warm (28°C) ambi- 
ent temperature (RHEN & LANG 1999). In addition, lizards 
mature earlier during longer hot seasons (ADOLPH & POR- 
TER 1996). Other factors should also be considered, such 


as the availability of food resources. Their contribution to 
the individual growth rate in F. labordi should be quan- 
tified to separate direct effects of temperature on growth. 
The current study is also limited by the fact that we did not 
follow individuals over time, but rather considered popu- 
lation averages. 

In conclusion, the current study suggests that climatic 
variation induces variation in hatching time and juvenile 
growth rate of chameleon species. A delayed rainy seasons 
induced by climate change promote late hatching, prompt- 
ing neonates to grow faster in compensation. Such chang- 
es probably affect embryonic diapause and hatching suc- 
cess. Annual species, such as F. labordi might be particu- 
larly vulnerable to climate change as a whole population 
could be lost if a single cohort cannot adjust their phenol- 
ogy with the year-to-year climatic fluctuations. Since the 
population in Ranobe has much shorter lifespan (KARSTEN 
et al. 2008), it is highly likely that southern populations 
will vanish first. Furcifer labordi is an extreme example for 
short lifespan among tetrapods, but other chameleon spe- 
cies seem to have an annual life history as well (ANDREONE 
et al. 2005, TOLLEY et al. 2010, GLAW et al. 2020). Therefore, 
several chameleon species, especially those from highly 
seasonal habitats in western Madagascar, may suffer from 
unexpected extreme weather conditions and might even 
become extinct due to temporary non-appearance of rainy 
seasons. 
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